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Abstract We present a variational method for deriving relativistic two-fermion wave equations 
in a Hamiltonian formulation of QED. A reformulation of QED is performed, in which covariant 
Green functions are used to solve for the electromagnetic field in terms of the fermion fields. The 
resulting modified Hamiltonian contains the photon propagator directly. The reformulation permits 
■ one to use a simple Fock-space variational trial state to derive relativistic fermion-antifermion wave 

equations from the corresponding quantum field theory. We verify that the energy eigenvalues 
' obtained from the wave equation agree with known results for positronium. 

in ' 

q ■ 1. Introduction 

m ; 

The description of relativistic bound and quasi-bound (i.e. unstable) few body systems continues 
to be an active area of research. The traditional method of treating relativistic bound states in 
quantum field theory (QFT) is by means of the Bethe-Salpeter (BS) equation. However, this 
approach has a number of difficulties, including the appearance of relative-time coordinates and 
• negative-energy solutions. In practice, the interaction kernels (potentials) in the BS equation are 

obtained from covariant perturbation theory, which may be of questionable validity for strongly 
coupled systems. In addition, the BS formalism is difficult to implement for systems of more than 
two particles. 

An alternative approach might be the variational method, which is non-perturbative in prin- 
ciple. The variational method has not been widely used in quantum field theory, in contrast to 
non-relativistic systems describable by the Schrodinger theory, in part because of the difficulty of 
constructing realistic yet tractable trial states. 

It has been pointed out in previous publications [1,2] that various models in QFT, including 
QED, can be reformulated, using mediating-field Green functions, into a form particularly conve- 
nient for variational calculations. This approach was applied recently to the study of relativistic 
two-body states in the scalar Yukawa (Wick-Cutkosky) theory [3,4,5]. In the present paper we shall 
implement this approach to the realistic QED theory, where comparison with experimentally ver- 
ified results are possible. In particular, we shall use the reformulated QED Hamiltonian to derive 
a relativistic fermion-antifermion wave equation and discuss its solution. 

The reformulation of the QED is presented in section 2, while the Hamiltonian and equal time 
quantization are given in section 3. In section 4 we use the variational principle, with simple Fock- 
space trial states to derive the relativistic fermion-antifermion equations, and present their "partial 
wave" decomposition for all possible J states. The relativistic radial equations are presented in 



X 
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section 5, while their non-relativistic and semi relativistic limits are given in 6. In section 7 the 
energy eigenvalues are shown to yield the correct the fine and hyperfine structure for all states. 
Concluding remarks are given in section 8. 



2. Reformulation of field equations and Lagrangian 

The Lagrangian of QED is (h = c = 1) 

£ = ~$(x) (i^d^ - m - e^A^x)) ip(x) 

~ {d a A p {x) - dpA a {x)) (d a A?{x) - d^A a {x)) . (1) 

The corresponding Euler-Lagrange equations of motion are the coupled Dirac-Maxwell equations, 

(iVfy - m) ^(x) = erfA^xMx), (2) 

and 

d^A"(x)-d»d^(x)=f(x), (3) 

where 

j"(x) = eip(x)~/ u ip(x). (4) 

The equations (2)-(3) can be decoupled in part by using the well-known formal solution [6,7] of the 
Maxwell equation (3), namely 

A^x) = A%x) + j d'x'D^x - x')f(x'), (5) 

where D^ u (x — x') is a Green function (or photon propagator in QFT terminology), defined by 

d a d a D^ u {x - x') - d^d a D au {x - x') = g^ v 5 A {x - x'), (6) 

and A°„{x) is a solution of the homogeneous (or "free field") equation (3) with j tl (x) = 0. 

We recall, in passing, that equation (6) does not define the covariant Green function D fJiU (x — x') 
uniquely. For one thing, one can always add a solution of the homogeneous equation (eq. (6) with 
9nv — ► 0). This allows for a certain freedom in the choice of D^, as is discussed in standard texts 
(e.g. ref. [6,7]). In practice, the solution of eq. (6), like that of eq. (3), requires a choice of gauge. 
However, we do not need to specify one at this stage. 

Substitution of the formal solution (5) into equation (2) yields the "partly reduced" equations, 



{i^df, - m) ij){x) = (a^x) + j d 4 x'D^{x - x')f{x')^ ip(x) 



(7) 



which is a nonlinear Dirac equation. To our knowledge no exact (analytic or numeric) solution of 
equation (7) for classical fields have been reported in the literature. However, approximate solutions 
have been discussed by various authors, particularly Barut and his co-workers (see ref. [8,9] and 
citations therein). In any case, our interest here is in the quantized field theory. 

The partially reduced equation (7) is derivable from the stationary action principle 

SS [il>}=5 j d A xC R = (8) 
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with the Lagrangian density 

C R = V^(x) {i-fd^ - m - e 7 ^(x)) Mx) - \ J d A x'f (x')D^(x - x')f{x) (9) 

provided that the Green function is symmetric in the sense that 

D^ u (x - x') = D^ v (x' - x) and D^ u (x - x') = D u ^(x - x'). (10) 

One can proceed to do conventional covariant perturbation theory using the reformulated QED 
Lagrangian (9). The interaction part of (9) has a somewhat modified structure from that of 
the usual formulation of QED. Thus, there are two interaction terms. The last term of (9) is 
a "current-current" interaction which contains the photon propagator sandwiched between the 
fermionic currents. As such, it corresponds to Feynman diagrams without external photon lines. 
The term containing Aq corresponds to diagrams that cannot be generated by the term containing 
D^ u , including diagrams involving external photon lines (care would have to be taken not to double 
count physical effects). However, we shall not pursue covariant perturbation theory in this work. 
Rather, we shall consider a variational approach that allows one to derive relativistic few-fermion 
equations, and to study their bound and scattering solutions. 



3. Hamiltonian of the quantized theory in the equal-time formalism 

We consider this theory in the quantized, equal-time formalism. To this end we write down 
the Hamiltonian density corresponding to the Lagrangian (9), with the term for the free Aq(x) 
field suppressed since it will not contribute to the results presented in this paper. The relevant 
expression is: 

Hr = Hq+Hi, (11) 

where 



Ho = ^{x){-i~a ■ V + mf3)^{x), 
Hi = ^Jd A x'f(x')D^(x-x')f(x), 



(12) 
(13) 



We construct a quantized theory by the imposition of anticommutation rules for the fermion 
fields, namely 

(x-y), (14) 



{v a (x,*),^(y»*)} =<W<5 3 



while all other vanish. In addition, if Aq ^ 0, there would be the usual commutation rules for the 
Aq field, and commutation of the Aq field operators with the field operators. 

To specify our notation, we quote the usual Fourier decomposition of the field operators, namely 



y,(*) = ]T 



d 3 p 



(27T 



,3/2 



rn 



to. 



1/2 



b ps u(p,s)e- i P- x + d\ )S v (p, a )e*- 



(15) 
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with p = p^ = (ujp, p), and uj p = y m 2 + p 2 . Dirac spinors u and u for free particles of mass m, 
where ("f^p^ — m)u (p, s) = 0, (7^^ + m) i> (p, s) = 0, are normalized such that 

u f (p, s) u (p, <r) = t» t (p, s)u(p,(7) = ^(5 Sffl (16) 
n t (p,s)w(p,a) = u f (p,s)u(p,cr) = 0. (17) 

The creation and annihilation operators 6' , b of the (free) fermions of mass m, and cr, d for the 
corresponding antifermions, satisfy the usual anticommutation relations. The non-vanishing ones 
are 

{b P s, b\ a ) = [d ps , d\ a ] = 6 sa 6 3 (p - q) . (18) 



4. Variational principle and fermion-antifermion trial states 

Unfortunately we do not know how to obtain exact eigenstates of the Hamiltonian (11). There- 
fore we shall resort to a variational approximation, based on the variational principle 

6(^\H-E\^ t=0 = 0. (19) 

For a fermion-antifermion system , the simplest Fock-space trial state that can be written down in 
the rest frame is 

IV-t) = E / d 3 pF slS2 (p)bl Sl dl pS2 |0> , (20) 

where F SlS2 are four adjustable functions. We use this trial state to evaluate the matrix elements 
needed to implement the variational principle (19), namely 

<V r | :H -E: ty T ) = £ j d 3 pF; iS2 (p)F SlS2 (p) (2u p - E) (21) 

S1S 2 J 

and 

( 27F ) s[s' 2 s lS2 J W ^ 

f -u{p,si)^u(p l : s' l )D llu {p-p l )v(-p l ,s l 2 )Yv(-P-,S2) \ (t)l) , 
X V +u(p,s 1 )^v(-p,s 2 )D^(p + p')v(-p',s' 2 )Yu(p',s' 1 ) J ' 

where p = (u> p , p), p' = (ay , p'), with p + p' = (i.e. p + p' = (2co p , 0)) in the rest frame, and 

D,u(x - x>) = j ^D^{k)e- ik ^-*'\ (23) 



We have normal-order the entire Hamiltonian, since this circumvents the need for mass renor- 
malization which would otherwise arise. Not that there is difficulty with handling mass renormal- 
ization in the present formalism (as shown in various earlier papers; see, for example, ref. [10] and 
citations therein). It is simply that we are not interested in mass renormalization here, since it 
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has no effect on the two-body bound state energies that we obtain in this paper. Furthermore, the 
approximate trial state (20), which we use in this work, is incapable of sampling loop effects. Thus, 
the normal-ordering of the entire Hamiltonian does not "sweep under the carpet" loop effects, since 
none arise at the present level of approximation, that is with the trial state | Y> T ) specified in eq. 
(20) . 

The variational principle (19) leads to the following equation 
J d 3 p(2u Jp -E)F SlS2 (v)5F; iS2 (v) 

, xs E /^ S ^ 1 «(q)H)-M.i««n«(P.q)^ 1 «(P) = 0, (24) 



S\S2 

m 



where M. SlS2<Tlcr2 (p,q) is an invariant "matrix element", which contains two terms: 

M SlS2aia2 (p, q) = M°ll 2(JirJ2 ( P , q) + M™ aiaa (p, q) , (25) 

where 

MZ e S2ai(T2 (p, q) = -u (p, si) (~ie^) u (q, o\ ) iD^p - q)v (-q, u 2 ) {-ierf) v (-p, s 2 ) , (26) 



>CW 2 (P, q) = « (P, si) (-icy") u ("P, «2) *iV (p + g) U (-q, <r 2 ) (-iej v ) u (q, <7i) , (27) 

correspond to the usual one-photon exchange and virtual annihilation Feynman diagrams. 

At this point it is worthwhile to make a few comments about our equation (24) and to compare 
its general features with other two-fermion equations, particularly field-theory based approaches. 
Firstly we note that the present variational derivation leads to momentum-space Salpeter-like 
equations, with at most four independent components F S1S2 (p). The equations have only positive- 
energy solutions, as is evident from eq. (24) with the interaction turned off, in which case only 
E = 2uj p > is obtained. This is in contrast to the BS equation, which is a 16-component equation 
and contains both positive, negative and mixed energy solutions. 

The interaction kernels, represented by the covariant A^-matrices, result from the variational 
derivation, that is, they are not put in by hands. This is in contrast to two fermion equations, which 
are not derived from a underlying quantum field theory, such as various two-body generalizations 
of the one-body Dirac equation. There are many such equations on the market, for example the 
eight component two-fermion equation of Pilkuhn [11]. In these treatments QFT effects, such as 
the virtual annihilation interaction (eq. (27)) do not arise naturally but need to be added in. 

The fact that only the lowest order ("tree level") diagrams appear in our eq. (24) is a reflection 
of the fact that we have used the simplest possible variational ansatz (20). Even so, it is important 
to note that, because of the reformulation discussed in section 2 and 3, their derivation does not 
require additional Fock-space terms in the variational state (20) as is the case in traditional (non- 
reformulated) treatments (e.g. [12]- [14]). 

In the non-relativistic limit, the functions F SlS2 can be written as 

F SlS2 (p) = F(p)A SlS2 , (28) 
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where the non-zero elements of A.^ for total spin singlet (S = 0) states are A12 = — A21 = A=, while 
for the spin triplet (S = 1) states the non-zero elements are An = 1 for m s = +1, A12 = A21 = 
for m s = 0, and A22 = 1 for m s = —1. We use the notation that the subscripts 1 and 2 of A 
correspond to m s = 1/2 and m s = — 1/2 (or f and j) respectively. Substituting (28) into (24), the 
variational procedure, after multiplying the result by A S1S2 and summing over s\ and s 2 , gives the 
equation 

(2u p - E)F(p) = j^j A^p, q)F(q), (29) 

where 



/C(p,q) 



m 



(p, q) A 



(71(72 • 



(30) 



S1S2CT1CT2 



To lowest-order in |p|/m (i.e. in the non-relativistic limit), the kernel (30) reduces to K, 
e 2 / |p — q| 2 , and so (29) reduces to the (momentum-space) Schrodinger equation 



2fi 



Hp) 



(2k) 



d 3 . 



Ip - q' 



(31) 



where e = E — 2m and fi = m/2. This verifies that the relativistic two-fermion equation (24) has 
the expected non-relativistic limit. 

In the relativistic case we do not complete the variational procedure in (24) at this stage to 
get equations for the four adjustable functions F SIS2 , because they are not independent in general. 
Indeed we require that the trial state be an eigenstate of the total angular momentum operator (in 
relativistic form), its projection, parity and charge conjugation, namely that 



(32) 



where mj = J, J— 1, — J as usual. We present explicit forms for the operators J 2 , J3 in Appendix 
A. The form for J 2 , eq. (109), in particular, is not readily found in standard texts and reference 
books. 

The functions F SlS2 (p) can be written in the general form 



" J 2 " 




" j(j + i) " 




V 


IY>T> = 


mj 
P 


\ij) T 


c 




C 





f SIS2 (p) = J2 E /£ 



S2 



(33) 



where yY 1 " 1 " 2 (p) are the usual spherical harmonics. Here and henceforth we will use the notation 

1 s 2 

p = |p| etc. (four-vectors will be written as p^). The orbital indexes £ slS2 and m slS2 depend on the 

spin indexes s\ and $2 and are specified by equations (32). The radial coefficients /si^ 2 ™ 8182 (p) m 
the expansion (33) also depend on the spin variables. 

Substitution of (33) into (20) and then into (32) leads to two categories of relations among 
the adjustable functions, as shown in Appendices A and B. It follows that, for trial states of the 
form (20), the total spin of the system is a good quantum number, and the states of the system 
separate into singlet states with the total spin S = (parastates) and into triplet states with 
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S = 1 (orthostates). We should point out that this phenomenon is characteristic of the fermion 
antifermion systems, which are charge conjugation eigenstates, and does not arise for systems like 
/i + e~. 

The singlet states 

In this case £ SlS2 = £ = J, mu = m,22 = and m.12 = m 2 i = mj. The nonzero components of 
F slS2 (p) are F n (p) = F 12 (p), F n (p) = F 21 (p) and have the form 

F slS2 (p) = f^ J (p)Y^(p), (34) 

where the relations between f[^ J {p) and f^^ip) involve the Clebsch-Gordan (C-G) coefficients 
C^"^ 2 , that is 

f^ J (p)=C { J f j Jms ^f J (p), (35) 

as it shown in Appendix A. We see that the spin and radial variables separate for the singlet states 
in the sense that the factors fsts 2 J {p) have a common radial function f J (p). Thus, for the singlet 
states we obtain 

*WP) = C ( J f j Jms ^f J (p)Y^(p). (36) 

The C-G coefficients C^/™^ 2 have a simple form: C^) Jmi1 = C ( ^] Jm22 = 0, C^/™ 12 = 

~C < jmj Jm21 = 1 ( see Appendix A). Therefore for the singlet states we can write expression (20) in 
the explicit form 

Wt) = j d 3 pf J ( P )Y™'(p) (b^dl pl - 6t/_ pT ) |0> . (37) 

These states are characterized by the quantum numbers J,mj parity P = (— 1) J+1 and charge 
conjugation C = (— 1) J . As we can see, the quantum numbers £ (orbital angular momentum), and 
total spin S are good quantum numbers for the singlet states as well. The spectroscopical notation 
is 

The triplet states 

The solution of the system (32) for S = 1 leads to two cases (Appendix A), namely £ SlS2 = £ = J, 
for which 

^ 1S2 (p)=/iri j (p)^r is2 (p)> (38) 

and £ S1S2 = £ = J =F 1, for which 
where 

mu = mj - 1, mi2 = m 2 i=mj, m 2 2 = mj + 1. (40) 

The expressions for /f lS2 (f>) in both cases involve the C-G coefficients Cj^ ms for S = 1 listed in 
Appendix A, that is 
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where the index m s is defined as 

m s = +1, when m slS2 = m±i, 

m s = 0, when m SlS2 = rnn = rri2i, (42) 

m s = — 1, when m slS2 = nii2- 

Thus, for the triplet states with I = J 

F S1S2 (P) = cfJ^f J (p)Y^(p). (43) 

These functions correspond to states, which can be characterized by the quantum numbers J,mj, 
parity P = (— 1)^ +1 and charge conjugation C = (— 1) J+1 . The orbital angular momentum £, as 
well as the total spin S = 1, are good quantum numbers in this case. The spectroscopic notation 
for these states is 3 Jj. 

For the triplet states with t = J =F 1 we obtain the result 

wp) = ^2J J " 1)ms / J -Hp)C7 S2 (p)+^ +1)m v J+ Hp)^7;r 2 (p), (44) 

which involves two radial functions f J ~ l {p) and f J+l (p) corresponding to I = J — 1 and £ = J + 1. 
This means that I is not a good quantum number. Such states are characterized by quantum 
numbers J, mj, P = (— 1) J , charge conjugation C = ( — 1) J and spin S = 1. In spectroscopic 
notation, these states are a mixture of 3 (J — 1) 7 and 3 ( J + 1) j states. 

The requirement that the states be charge conjugation eigenstates (the last equation of (32)) is 
intimately tied to the conservation of total spin. Indeed, a linear combination of singlet and triplet 

F S1S2 (p) = Ci/iff? \p)Y^ (p) + C7 2 /iS J (p)Y^ (p), (45) 

satisfies the first three equations of (32). However, it is unacceptable for describing a fermion- 
antifermion system because the first and the second terms in (45) have different charge conjugation. 
For a system of two particles of different mass (such as u + e~) charge conjugation is not applicable, 
so that the total spin would not be conserved. 

5. The relativistic radial equations and application to positronium-like systems 

We return to equation (24) and replace the functions F SlS2 (p) by the expression (36) for singlet 
states and by (43) and (44) for triplet states. The variational procedure then leads to the following 
results: 

For the singlet states I = J, P = (— 1) J+1 , C = (— 1) J , the radial equations are 

(2^ - E) f J (p) = -f~ 3 [ ^/C(^) (p, q) f(q), (46) 

where the kernel 

JC^(p,q)= [dpd<iC^ s ^M SlS ^ (47) 
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is defined by the invariant M-matrix and the coefficients 

r i(sgl)siS2(Ti<T2 _ ^(sgtyJma n (sgl)Jm s , ST^ ( n {sgl)Jm v \ 2 ,,„•. 
U Jmj = U Jmj U Jmj I 2^ \ U Jmj ) " ^> 

u 1 i/ 2 m J 

Here we have summed over mj, because of the (2 J + l)-fold energy degeneracy. 

For the triplet states, we obtain different equations for the I = J, and I = J =F 1 cases. Thus 
for the states with t = J, P = (-1) J+1 , C = (-1) J+1 the result is 

(2u p - E) f\p) = -^3 / ^JC^\ P , q )f J (q), (49) 
where the kernel if( tr ) is formally like that of (47), namely, 



/C< fr >(p,g)= C { jZ lS2aia2 j dvd^M^Av^YT^iWT 1 * 2 (q). 



(50) 



However it involves different C-G coefficients, namely 



C 



(tr)siS2<T\(T2 

Jmj 



(51) 



v\V2mj 



For the triplet states with I = J^fl, we have two independent radial functions f J ~ l {p) and f J+1 {p)- 
Thus the variational equation (24) leads to a system of coupled equations for f J ~ 1 (p) and f J+l (p)- 
It is convenient to write them in matrix form, 



where 



and 



(2co p - E) F (p) = JpL / £!*K(p >g )F(g), 

F(p): 



/Cn (p, g) /C12 (p, g) 
(p, g) £ 22 (p, q) 



K(p,q) = 

The kernels Kij are similar in form to (47) and (50), that is 



(52) 
(53) 
(54) 

(55) 



IC ij (p,q)= CjZT 2 /dpcZqM SlS2CTlCT2 (p,q)y,7 CT1CT2 (q)^7 slS2 *(p). 

CTl(T2SlS2mj 

However the coefficients C S j^^ U2 are defined by expression 

Jmjij Jmj Jm.j ' / j \ Jm.j J ' 

siS2mj 

where l\ = J — 1, £2 = J + 1 and ms is as defined in Eq. (42). The system (52) reduses to a single 
equation for J = since / J_1 (p) = in that case. 



(56) 
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Our equation (24), or its radial components (46), (49), (52), contain the relativistic two-body 
kinematics (kinetic energy, recoil effects) exactly, but the dynamics are included approximately 
due to the limited nature of our trial state (20). This limitation is reflected in the fact that the 
interaction kernels of our equations contain only "tree-level" Feynman diagrams. Nevertheless our 
equations (46), (49), (52) have no negative-energy solutions, in contrast to the BS equation. They 
are variationally derived, hence the energy eigenvalues obtained from them will give meaningful 
values for any strength of the coupling. 

To our knowledge, it is not possible to obtain analytic solutions of the relativistic radial 
momentum-space equations (46), (49) and (52). Thus one must resort to numerical or other 
approximation methods. Numerical solutions of such equations are discussed, for example, in [10], 
while a variational approximation has been employed in [5]. However, in this paper we will concen- 
trate on perturbative 0(q 4 ) solutions, since it is important to verify that our equations yield the 
correct fine structure for systems like positronium. 

Our equations will yield energies which are incomplete beyond 0(a 4 ), because our variational 
trial state (20), as mentioned, reflects only "tree- level" Feynman diagrams, that is no radiative 
corrections are incorporated. One could, of course, augment them by the addition of invariant 
matrix elements corresponding to higher-order Feynman diagrams (including radiative corrections) 
to the existing .M-matrices in the kernels of our equations, as is done in the BS formalism. Indeed, 
such an approach has been used in a similar, though not variational, treatment of positronium 
by Zhang and Koniuk [15]. These authors show that the inclusion of invariant matrix elements 
corresponding to single-loop diagrams yields positronium energy eigenstates which are accurate to 
O (a 5 , a 5 In a). However such augmentation of the kernels "by hand" would be contrary to the 
spirit of the present variational treatment, and we shall not pursue it in this work. 



6. Semi-relativistic expansions and the non-relativistic limit 



For perturbative solutions of our radial equations, it is necessary to work out expansions of 
the relevant expressions to first order beyond the non-relativistic limit. This shall be summarized 
in the present section. We perform the calculation in the Coulomb gauge, in which the photon 
propagator has the form [16] 



A,o(k) = ^, Ay(k) = 0, D ij (k») = - j J^ 



ki kj 

on - -gr 



(57) 



where fc M = (u p — uj q , p — q). 

To expand the amplitudes Ai of (26) and (27) to one order of (p/m) 2 beyond the non-relativistic 
limit, we take the free-particle spinors to be 



u(p,i) 



1 + ^ 

C?-P) 
2m 



<Pi, 



v(p,i) 



2m 
x ~ Srn? 



(58) 



as discussed in Appendix C. In this approximation the photon propagator takes on the form 

(p ~ q)i ip - q] 



A)0 (P 



(p-q) 



■ — a> Aj(p-q)-— 



(p - q)' 



Sij - 



(p - qr 



(59) 
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Corresponding calculations give for the orbital part of the .M-matrix 

M ope{orb) , x = - e 2 f 1 J_ A _Q_P_ (p_Xq£\ 1 , . 

The terms of the expansion linear in spin correspond to the spin-orbit interaction: 
M o P e(s-o) , } = j£» t t ~ • (P X q) 

ip — qj 

Here "ir ( + ) and "a ( _ ) are positron and electron spin matrices respectively, defined as follow: Tj ( + ) (p ai x s 

( 0^ + V<ti) X S2 > "<^ < '~Vo-iXs2 = Wi (^~^Xs 2 )- ^he Quadratic spin terms or spin-spin interaction 
terms are 



;„2 



ie 



vUU < — v " " /yv ; 2 — + i?w . i ^ iXs2 . (62) 



(^ (+) • (P - q)) • (P - q)) . -,(_) 



4m 2^lA<X2| (p _ q) 

Lastly, the virtual annihilation contribution is given by 

M™ nn {p t q) = -^Ul 2 {^ (+) • ^X S2 , (63) 

where we have excluded a divergent term, which appears in the Coulomb gauge calculation. This 
divergence is an artefact of the Coulomb gauge. It does not arise, for example, in the Lorentz 
gauge, where only the expression (63) is obtained. However the Lorentz gauge is not convenient for 
obtaining all other 0(a 4 ) corrections because it contains spurious degrees of freedom (longitudinal 
polarization) of the photon. 

We have used expressions (60)-(63) to obtain the corresponding radial kernels. Details of the 
calculations can be found in Appendix D. We use the notation z = (p 2 + g 2 ) /2pq, and Q\(z) is 
the Legendre function of the second kind [17]. The contributions of the various terms to the kernel 
are as follows: 

Singlet states with £ = J (J > 0), P = (-1) J+1 , C = (-1) J 
Orbital term 

)C(s9l)(o) (p >g ) 

2we2 Qj(z) + ^4 (~ + ~) QM + ( J + i) Qj+M ~ 2^,ol • (64) 



pq m 2 \ 2 \q p 

Spin-orbit interaction 



Spin-spin interaction 



K,^ s -°\p,q)=0. (65) 



K {sgt){s-s)^ ) = 2EL. 8 (66) 
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Triplet states with £ = J (J > 1), P = (-1) J+1 , C = (-1) J+1 
Orbital term 

K^Hp, q) = —Qj (z) + ^ (~ + -) QJ (*) + ^ + 1) (*)) • (67) 

m z \ 2 \q p J J 



Spin-orbit interaction 

K (tr)(s-o) {pj g) = {Qj+l {z) _ Qj l {z)} . (6g) 



Spin-spin interaction 

= £ (f + 1) Q > < 3 » - S 27TT ( JQj +' <*> + < J + » <- )} ' (69> 

Triplet states with £ = J - 1 (J > 1), £ = J + 1 (J > 0), P = (-1) J , C = (-1) J 

The off-diagonal elements of the kernel matrix (Eqs.. (52)- (54)), K\2 and K21 which are 
responsible for mixing of states with £ = J — 1 and £ = J + 1, get a non-zero contribution from the 
spin-spin interactions only: 

/C 12 (p, g) = JC 21 (p, q) = (jQj+i (z) + \Qj-i (z) - 2Qj (*)) . (70) 

The contributions to the diagonal elements of the kernel matrix are the following: 
Orbital terms 

lC$(P,q) = —Qj-i(z) 
pq 

+^ (-^ (f + 1) Qj-i (*) + JQj (*) - 2*7-1,0) , (7i) 



m z \ 2 \ q p / / 

Spin-orbit interaction 

*fr o) (p, ?) = wj oo - qj-2 w) , (73) 

4r° } M = - ^ 2 j + 3 (z) - Qj (*)) • (74) 
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<T S) (P, ?) = S^ttt^t ((z + z) Qj-i (*) ~ 2 Qj (*) ) > ( 75 ) 

42" S) (P,5) = ^4:77^ (Y? + I ) (*) " 2Qj + 2 (z) ) • (7(5) 



Spin-spin interaction 

\ 

in /Tl — - - - 

2m 2 2J + 1 \\q p 

ve 2 1 ffP.g 
2m 2 2J + 3\\q p 

Annihilation term 

2ire 2 

JC ann ( P ,q) = r 8j- lfl . (77) 

We note that in the non-relativistic limit only the first terms of the orbital part of the kernels 
survive. They have the common form 2irie 2 Qe(z)/pq, hence all radial equations reduce to the form 

m 2 2 roc 

(2oj p - E) f(p) = / dq±Q e (z)f e (q). (78) 

TTUJpP J UJ q 

Recalling, also, that 

Wp = v/m^ + p 2 ~ m { 1 + l(Z.y\ (79) 
we obtain, in the non relativistic limit, the momentum-space Schrodinger radial equations 

(^-e)f\p) = -- f°° dqqQj(z)f e (q), (80) 

where a = e 2 /Air, /x = ^, e = E — 2m. 

7. Energy levels. Fine and hyperfine structure 

The relativistic energy eigenvalues E n j can be calculated from the expression 

roc roo 

E dpp 2 f J (p)f J (p) = / dpp 2 2io p f J (p)f J (p) 
Jo Jo 

m 2 t°° dpp 2 r °° 



(2vr) d Jo w. 



— r dq^-K,^ tr \p,q)f J (p)f J {q) (81) 

J p Jo u q 



for the singlet and I = J triplet states. 

For the t = J T 1 triplet states the corresponding result is (see equation (52)) 

POO poo 

E / dpp 2 ¥ jf (p)¥(p) = / dpp 2 2uj p ¥\p)¥{p) 
Jo Jo 



m r dpp_ f°° dq ^ KM¥J{pnq y (g2) 

2-K) Jo Up Jo LO q 



To obtain results for E to O (a 4 ) we use the forms of the kernels expanded to O [p 2 /m ) (eqs. 
(64)-(77)) and replace / (p) by their non-relativistic (Schrodinger) form (see (147), Appendix D). 
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The most important integrals that we used for calculating (81) and (82), are given in Appendix 
D. In Appendix E we show that the contribution of kernels if 12 and K21 in (82), is zero at 0(a 4 ). 
Thus, the energy corrections for the triplet states with £ = J — 1 and £ = J + 1 can be calculated 
independently. 

The results will be presented in the form Ae = E — 2m + a 2 m/2n 2 . 

Singlet states (l = J (J > 0), P = (-1) J+1 , C = (-1) J ) 
The kinetic energy corrections 

Ae^ = -^(-?— V !J zV (83) 
K 8 \2J + ln 3 8n 4 J K 1 

The potential energy corrections 



Xsgl){o) ® 4 m ( ( 3 ^ \ 1 I 



=-^ll2JTT- 2 H^-^ • < 84 > 



A£ (sgl)(s-o) = Q) (85) 



Ae ^-) = -^% (86) 
^ 4 n 3 



The total energy corrections 



8 \2J + ln 3 8n 4 y 



Triplet states (i = J (J > 1), P = (-1) J+1 , C = (-1) J+1 ) 
The kinetic energy corrections 



,., a 4 m ( 1 i 3 1, 

*' = 8~~ ( 2jTTi? ~ 8 n* ) ' <88> 

The potential energy corrections 



Ae ^> = (9i) 



_( s _ s ) _ arm 

8 J(J + l)(2J + l)n 3 ' 

The total energy corrections 



(87) 



^"•»=-=HG^-**)"], (89> 
A4 ,,-^_^ _j»_ i, (90) 



Ae (-) = (Y_L_ + ? ^ J_ _ 11 J_ ) (92) 

8 VW + l «/(</ + l)(2J + l)/n3 8n 4 /' 
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Triplet states (l = J - 1 (J > 1), P = (-1) J , C = (-1) J 
The kinetic energy corrections 



A (tr)(J-i) _ _ 



4 

am 



1 1 3 1 



2 J - 1 n 3 8 n 4 



The potential energy corrections 



Ae 



(tr)(o)(J-l) 



4 

am 



2J- 1 



2<5j,i 



1 1 " 

n 3 n 4 



a tr(s-o)(J-i) a 4 m3(l -8 Jt i) 1 
p 8 J(2J-l)n 3 ' 



Ae 



tr(s-s)(J-i) _ am 



1 — 6 



J A 



8 J (2J + 1) (2J - l)n 3 ' 



4 n 3 ' 



The total energy corrections 



Ae tr(J-l) 



4 

am 



4 _ 2 (3J + 1) (1 - fyi) 
2J-1 J(2J + 1) (2J- 1) 



Triplet states ^ = J + 1 (J > 0), P = (-1) J , C = (-1) J ) 
The kinetic energy corrections 



A (tr)(J+i) = _<a_m / 1 1 _ 3 J_ 
K 8 \2J + 3n 3 8n 4 



The potential energy corrections 



Ae (tr)(o)(J+l) = 



m 



3 1 1 



2 J + 3 n 3 n 4 

a tr(s-o)(J+i) _ « 4 m 3 1_ 

£p 8 (J + l)(2J + 3)n 3 ' 



£p 8 (J + l)(2J + 3)(2J + l)ra 3 ' 



The total energy corrections 



Ae 



tr(J+l) _ ^ m ( 2 • ! - / + 2 



1 11 



2J + 3 V (J + 1) (2J+ 1)7 n 3 
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These results are in agreement with the well-known positronium fine structure results [18], [19]. 



8. Concluding remarks 

We have considered a reformulation of electrodynamics, in which covariant Green functions are 
used to solve the field equations for the mediating electromagnetic field in terms of the fermion field. 
This leads to a reformulated Hamiltonian with an interaction term in which the photon propagator 
appears sandwiched between fermionic currents. 

The variational method within a Hamiltonian formalism of quantum field theory is used to 
determine approximate eigensolutions for bound relativistic fermion-antifermion states. The refor- 
mulation enables us to use the simplest possible trial state to derive a relativistic momentum-space 
Salpeter-like equation for a positronium-like system. The invariant M matrices corresponding to 
one-photon exchange and virtual annihilation Feynman diagrams arise directly in the interaction 
kernel of this equation. 

The trial states are chosen to be eigenstates of the total angular momentum operator J 2 and J3, 
along with parity and charge conjugation. A general relativistic reduction of the wave equations 
to radial form is given. For given J there is a single radial equation for total spin zero singlet 
states, but for spin triplet states there are, in general two coupled equations. We show how the 
classification of states follows naturally from the system of eigenvalue equations obtained with our 
trial state. 

It is not possible, as far as we know, to obtain analytic solutions of our relativistic radial 
equations nor the resulting eigenvalues of the particle-antiparticle system described. However, it 
is possible to obtain 0(a 4 ) corrections analytically for all states using perturbation theory. The 
results agree with well known results for positronium, obtained on the basis of the Bethe-Salpeter 
equation [19], which lends credence to the validity of our variationally derived equations. 

The method presented here can be generalized to include effects higher order in alpha by using 
dressed propagators in place of the bare propagators. This shall be the subject of a forthcoming 
work. 
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Appendix A. Total angular momentum operator in relativistic form 

The total angular momentum operator is defined by expression 

J = J d 3 x:^ (x) (L + S)V> (x) :, (104) 

where L = x x p and S = }fa are the orbital angular momentum and spin operators. We use the 
standard representation for the Pauli matrices 



a 



a 
~a 



(105) 
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" 






" -i ' 




" 1 


1 





0-2 = 


i 


, 0-3 = 


-1 



(106) 



. we obtain 




,t f, . i At 



L q2 (&q T &qT + 6 q| & Ql + 4?^ + 4^) 
I + & q| ft qT " 4t^i + 4^, 

Lq3 (& qT &qT + & q| 6 ql + 4t^t + 4l^ 



(107) 



% +5 ( & qT fe qT - 6 q| 6 ql + 4?^ " d q^, 

Here L g is the orbital angular momentum operator in momentum representation: 

(L,) i = L 9i = -i(qxV,) i . (108) 

Note that these expressions are valid for any t, since the time-dependent phase factors of the form 
e iu> q t cance i ou t. 

For the operator J 2 = j\ + j| + ^3 we have 



J 2 = 



v I + I ) ( & W + & q|&ql + 4^ + ^ql) ^ 



+ 



/ 



2L q t ■ L 



& q'T Vt^qT^T + 6 q' T Vt^q^qt 

- y +^ i Vi d qT d qT + fe q'|Vi d qi d qt 
+\ ( &J, t Vt4t^t - 4'tVt4i^i 

+ 6 ^t Vi4i^t + 4i Vt4t^ 

4t Vr4t^t + 4i Vi4i d qt 
4^4'^ + 6 Vqt4'A'i 
4t Vt4t d it + 6 q'i Vi4T d qt 
+4 T Vi4A't + 6 qt 6 qi4'A'i , 

^ q , T 6 q / T 4 T ^qT - frq'l Vt4t^t 




where 



£<2+ = ^gl + iL q 2, L q - = L q \ — lL q 2 , 



'qt 



(109) 



(110) 
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The requirement that the trial state (20) be an eigenstate of J 2 and J z leads to the system of 
equations 

(l 3 + i)Fh = mjFu, 
L3F12 = mjFu, 

L3F21 = mjF 21 , (111) 
(Z 3 - l) F 22 = mjF 2 2, 

( J( J + 1) - L 2 - 2 - 2L 3 ) F u = L_(F 12 + F 21 ), 

( J( J + 1) - L 2 - l) F12 = F 21 +L + F 11 +L-F 2 2, 
( J( J + 1) - L 2 - l) F 2 i = Fis + L+Fn + L.^, (112) 
( J( J + 1) - L 2 - 2 + 2L 3 ) F 22 = L+(F 12 + F 21 ). 
Substitution of the expressions (33) for F S1S2 and use of eq. (Ill) gives 

mu = m2i = rnj, run = mj — 1, 11122 =mj + l, (113) 



'-11 = £22 



£12 = £21 = e, (ii4) 



and 



(j(j + i)-e(e + i)-2m J )fl 1 (p) = V( e-mj + i)(e + ^j)fi 2 {p) 

+V(£-mj + l)(£ + mj)f% 1 (p), 
{ j {J + 1) _ e{e + 1) _ 1)f t 2{p) = f 2i(p ) 

+ V(£ + mj)(£-mj + l) f n (p) 

+ V(£-rnj)(£ + mj + l)fZ 2 (p) , 

{j(j + i)-e(e + i)-i)f^(p) = f{ 2 (p) (115) 

+ V(£ + mj)(£-mj + l) f e n (p) 

+ y(e-mj)(e + m J + i )fi 2 (p) , 

(J(J+l)-e(e + l)+2mj)fZ 2 (p) = y/(t + rnj + l)(e-mj)fi 2 (p) 

W(£ + mj + m-mj)f 2 i{p). 

The singlet states correspond to the solution f[ x (p) = / 22 (p) = 0, ff 2 (p) = —f£i (p) °f this 
system with £ = J (J > 0). 

For the triplet states the solutions are /f 2 (p) = / 21 (p) = f (p), and, 
for £= J- 1 (J> 1): 



{J-mj)fk\p) = y/(J -mj)(J + mj-l)f J -\p), (116) 
(J + mj)/^- 1 ^) = y/(J + mj)(J-m J -l)f J - 1 (p), (117) 
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for £ = J (J > 1): 



mj/ 2 J 2 (ri = VfJ-mjJfJ + mj + l)/^), 



for £ = J + 1 (J > 0): 

(J + l + mj )fJ+\p) 



■ V(J-mj + 2) (J + mj + l)/^+i (p) , 
- V(J-mj + l)(J + m J + 2)/ J+1 (p) . 



It is convenient to introduce the table of coefficients Cj^ m " : 



J -I 



£ = J 



£ = J+l 



TO, 



+ 1 



(J+mj-l)(J+mj) 
J(2J-1) 



(J+mj)(J-mj+l) 
J(J+1) 



(J-mj+l)(J-mj+2) 
(J+l)(2./+3) 



TO, = 



(J-mj)(J+mj) 
J(2J-1) 



(J-mj+l)(J+mj+l) 
(J+l)(2J+3) 



TO., 



(J-m;-l)(J-m/) 
J(2J-1) 



(J-mj)(J+mj+l) 
J(J+1) 



(J+mj+2)(J+m, 7 +l) 
(J+l)(2J+3) 



(118) 
(119) 



(120) 
(121) 



These coefficients coincide with the usual Clebsch-Gordan coefficients for S = 1 except for a 
factor 2 in the denominator, which we absorb into the normalization constant. 



Appendix B. Parity and charge conjugation 

We consider the application of the parity operator to the trial state (20): 

VWt) = E f d 3 P F SlS2 (p)Vbi Si dl pS2 \0) 

S1S 2 J 
S1S2 J 

Making use of the properties 

Vb\ Sl V- x = v P bl pSl , Vdl pS2 V- 1 = - V P 4 S2 , V |0> = |0) , 

where 7] P is the intrinsic parity ((r] P ) 2 = 1), it follows that 

V^t) = E/ d 3 P F SlS2 (p)Vbl Sl dl pS2 |0) = -E/ d 3 P F SlS2 (-p)bl Si dl pS2 |0) 



= PE/ d 3 pF SlS2 (p)bl Si dl pS2 \0) , 



(122) 



(123) 



(124) 



where the parity eigenvalue P depends on the symmetry of F SlS2 (p) in different states: 
For the singlet states (£ = J) we get from (36) F slS2 (-p) = (-1) J F SlS2 (p), so that P = (-1) J+1 . 
For the triplet states with £ = J we get from (38) F SlS2 (—p) = (— 1) J F slS2 (p), hence P = (—1) 



J+i 
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For the triplet states with t = J ± 1 we get from (39) F SlS2 (— p) = (—l) J+1 F SlS2 (p), therefore 
P=(-l) J . 

Charge conjugation is associated with the interchange of the particle and antiparticle. Applying 
the charge conjugation operator to the trial state (20) we get 



C|Vt) = E [ d 3 pF SlS2 (p)Cbl Sl dl pS2 \0) 

sis 2 J 
sis 2 J 

Using the relations 

a^c- 1 = r, c d\ sv Cdl pS2 c- 1 = n c bl pS2 , C |0> = |0) , 
where (r] C ^ 2 = 1, we obtain 

C Wt) = E / d 3 pF slS2 (p)Cbl Sl dl pS2 |0) = E / d 3 pF S2Sl ( P )bl Sl dl pS2 |0) 

Sl s 2 J 

= CE/ d 3 pF SlS2 (p)bl Sl dl pS2 \0) , 



(125) 
(126) 

(127) 



SlS 2 



(128) 



SlS 2 



where the charge conjugation quantum number C depends on the symmetry of F slS2 (p) in different 
states: 

For the singlet states (£ = J) we get from (36) F slS2 (—p) = (— 1) J+1 F slS2 (p), hence C = (— 1) J . 
For the triplet states with I = J we get from (38) F SlS2 (—p) = (— 1) 3 F SlS2 (p), therefore C = 

("1) J+1 . 

For the triplet states with £ = J ± 1 we get from (39) F SlS2 (—p) = (— 1) J+1 F SlS2 (p), so that 

c = { \y>. 



Appendix C. Expansion of the spinors 

We recall the form of the particle spinors: 

u(p,t) = N p 



1 

uip+m 



where 



1 




^2 





1 



ojp + m 
2m 



The antiparticle or "positron" representation for the Vi(p) spinors has the form 

v(p,i) = N p 



(?-p) 

LOp+m 
1 



where 



Xi 



X2 



Xi 
1 





(129) 
(130) 

(131) 
(132) 
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The normalization is 

u(p,i)u(p,j) = Sij, v(p,i)v(p,j) = -Sij. 
Expanding in powers of p/m and keeping the lowest non-trivial order terms, 

(o>-p) 



we obtain the result 



u(p,i) 
v(p,i) 



uj v + m 



2m 



~ 1 + 



~ r 8m 2 ) 
8m 2 J 



uj p + m ^ ^ p 



2m 



8m 2 ' 



1 

C?. P ) 

2m 




C?. P ) 

L 2m J 


2m 
1 




" C?-P) ] 
2m 

. > + & - 



(133) 

(134) 
(135) 

(136) 
(137) 



Appendix D. Some useful identities and integrals 

The following identity is useful for evaluating the Ai matrices: 



((p-q)-p) 



P 



(p x q) 



(p-q) 2 (p-q) 4 ' 



(P " q) 4 

The angular integration in (47), (50), (55) involves the following integrals 

J dpdHF (p • q) 3^' J (q)l7 J *(p) = 2ir6j,j6 m , jmj j d{p • q) F (p ■ q) Pj (p ■ q) , 



(138) 



(139) 



/ 



d (p ■ q) - P q . 2 ^J (p • q) 

(p-q) 



p| |q| \2J+ 1 



2J + 1 



(140) 



... mIpxqTp r ^ (J + l)(J + 2) ^ J(J-l) 

d(p " q) (p^ Pj (p " q) = 2( 2 j + i) {z) ~ WITT) 



Qj-i (z) , 



(141) 



where F (p • q) is an arbitrary function of p • q, Pj (x) is the Legendre polynomial, and Qj{z) is 
the Legendre function of the second kind of order J. 

The following integrals are needed for the calculation of the relativistic energy corrections. 



/*oo /*oo 2 

/ / dpd q pVf J (p)f J (q)=2n(^) Sj, 
Jo Jo \ n / 



(142) 
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poo poo 

/ / dpdqpqf J (p)f J (q)Qj(z 1 ) = ^, (143) 
Jo Jo n 

poo poo 

/ / dpdqp 2 q 2 f J (p)f J (q)Qj(z 1 ) = 
Jo Jo 

jfjf = »(v) 3 (wTT^)' < 144 > 

[° r dpd 4 p\ 2 f J ij>)f'{q)Qj- l (z 1 ) =tt (^)' 0- i~) , (145) 

jf" f*«i W ¥/ J (p)/'(« Jtl W = -r (f )* ( -J^ - I) . (146) 
Here J" 7 is the nonrelativistic hydrogen-like radial wave function in momentum space [19] 

J iP) - Jn iP) - ^ {n + J)l ) ^ + 1)J+2 y n -J-l [ n 2 p 2 + 1 )> ^) 

where G'^ l J _ 1 (x) are Gegenbauer functions. 

Appendix E. /C 12 , /C 2 i kernels for I = J =f 1 states 

The contribution of the kernel K\2 to the energy correction is 

J dpdqp 2 q 2 IC 12 (p,q) f J - 1 (p)f J+1 (q), (148) 

where 

/C 12 (p, g) = J] / dp dq^^) (p, q) (q)^ 2 * "(p) - (149) 

This requires the following integral 

E C jZ?2 U2 J ^pA/'" 1 ^^*^ (150) 

CT1CT2S1S2 

We calculate this form in coordinate space. The Fourier transform of M SlS2(7l(72 (p,q) is 

M SlS2ai(72 (p,q) = y d 3 rdV.M SlS2(Tl(T2 (r, r') e^P-iM'-^ (151) 
where the M SlS2(J1(J2 (r,r') matrix is a local operator in general [16], that is 

M SlS2IJl a 2 (r, r') = M SlS2aia2 (r) 5 (r - r') . (152) 
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We apply this transformation to the Msfs^aial (p ; q) matrix (see eq. (62)). Because of the angular 
integration in (149), only the first term in (62) survives. The Fourier transformation of that term 
is 

(*<+> ■ (p - q)) (*<-> ■ (p - q)) (gW ■ r) ( ?(-) ■ r ) 



Am 2 (p - q) 2 167rm 2 r 5 



Furthermore, 



j d\f J+ \ q )Y^\^)e-^ = R J n +1 (r)Y^r(r), (155) 



where 



The associated Laguerre function (p) is related to the confluent hyper geometric function by 

L " {p) = /x!( ( A-/i)! F ( " A + M ' ^ + 15 P) ' (157) 

The generating function for the Laguerre function is 

hence 

E / dWq/ J -W7T 2 *(p)^^ 

= E C s j^ J SvRi-\r)Y^\i) 



CT1CT2S1S2 



167rm 2 r 5 

3a 



167rm 2 

x 



| drr 2 !^- 1 ^)^ 1 (r)x 



(TltT2«lS2 

It follows that 

E <w$r I *Y£r\*) ( - (+) • • (*) = (i«>) 

fTl(T2SlS2 

but 



r°° i 

j( drr 2 -3^- 1 (r)^+ 1 (r)=0. (161) 
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The last expression can be proved in the following way. Let us consider the more general 

/»oo 

/ drr^Ri(r)RUr). 
Jo 

The generating function for R e n (r) is 



G ne (r,u)= £-e r/n — (-1) + ^xrrexp^-- 

n 2 Y ((n + ^)!) 3 VW (l-u) 2£+2 I 1-un 

Then we consider the expression 

/ drr^ +2 G^(r,u)G n ^(r, V ) 
J o 



'0 

coo 



4 \{n-t- l)!(n-f - 1)! /o_\<+^ 



ex P i - ( T— + 



(l-u) 2 * +2 (l-^' +2 



2r/n ^ 


1) 








\ 2r 


1-U 


/ n 



X 



2£+l„2£'+l 



_£ / (n-l- l)!(n-f -1)! u^+V _ 

n 4 V (( n + ^)!) 3 ((n + f)!) 3 (1 - +2 (1 - ^'+ 2 



exp ^ - 1 + + 



n J y \ 1 — u 1 — v J n 

It is well known that 

roo 

/ dppPe-P = r (/? + 1) , 
Jo 



therefore 



\ / w i> \ 2r 

exp { - 1 + + 



/ [ \ 1 — u 1 — v ) n 

p+3+e+e' 



and 



Jo V n 

poo 

/ drrP +2 G n£ (r,u)G n£ ,(r,v) 
Jo 



2-P- 1 / ( n -l- l)!(n-f - 1)! 
n-P +1 \l (( n + £)\) 3 ((n + £')\) 3 
u 2i + i v 2e> + i (1 _ u f+i-t+e' (1 _ 

X (i _ uv f+W 

We expand this expression in a series, 



T((3 + 3 + £ + £') . 



[■OO 

/ drr /3+2 G„, (r, u) G n& (r, v) = V C w (n, /3, £, (!) vPvP' 

Jo , 

vv 
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It is not difficult to show [20], that the coefficient C n+ £ iJl+ £i represents the integral 

poo 

C n+e , n +e(n,l3,£,t)= drr^ +2 R e n (r)R^(r). (169) 

J o 

Simple but tedious calculations show that this coefficient is zero for (3 = —3, £ = J — 1, £' = J + 1. 
Thus the kernel K\2 does not contribute to the energy corrections to O (a 4 ) . The same result is 
obtained for the kernel if 21- 
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